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We investigate three different methods to tackle the problem of diffusion-limited reactions (anni- 
hilation) of hard-core classical particles in one dimension. We first extend an approach devised by 
Lushnikov and calculate for a single species the asymptotic long-time and/or large distance behavior 
of the two-point correlation function. Based on a work by Grynberg et ai, which was developed 

Qto treat stochastic adsorption-desorption models, we provide in a second step the exact two-point 
correlation function (both for one and two-time) of Lushnikov's model. We then propose a new 
' formulation of the problem in terms of path integrals for pseudo-fermions. This formalism can be 

used to advantage in the multi-species case, specially when applying perturbative renormalization 
group techniques. 



PACS number(s): 05.70 Ln; 47.70-n; 82.20 Mj; 02.50-r 



I. INTRODUCTION 

C/3 , 

The recent interest in modelling low dimensional diffusion-limited reactions has been stimulated in part by the 
experimental observation of anomalous kinetics in low-dimensional systems Q . The traditional approach to chemical 
reactions is based on mean-field theory, i.e., rate equations for the densities of the various reactants. The latter 
I describe well the reaction kinetics in three dimensions because diffusive transport of reactants allows to eliminate the 
spatial fluctuations of the concentrations. However, in lower dimensions, due to the lack of phase space, the reactants 
spatial fluctuations can grow and develop inhomogeneities in the concentrations. Futhcrmore, even in the spatially 
homogeneous case, the rate equations are not applicable in less than three dimensions; for example, in the two-species 
diffusion-limited annihilation the concentration of the particles decays (for identical initial concentrations) slowlier 
than the mean-field theory predicts. Thus the fluctuation-dominated dynamics is beyond the classical theories, yet can 
^ ■ be accounted for by simple one-dimensional models of hard-core particles. The latter are solved either numerically 
or analytically by applying techniques from (classical or quantum) statistical mechanics In particular, exact 
solutions have been obtained by inter-particle distribution methods, by relating the systems to dual and solvable 
ID models [||J|-01 (Kinetic Ising and Potts) and by mapping the diffusion-reaction processes to an imaginary-time 
dynamics of quantum spin chains with non-hcrmitian Hamiltonians . An alternative and fruitful approach has been 
OO ' developed by Cardy and collaborators ||s|-|lTf : the idea is to reformulate the original problem in terms of a field theory 
Q\ , of interacting bosons and subsequently use renormalization group techniques. This is a powerful method as it applies 
■ to arbitrary dimension and low densities of particles, a regime where universal behavior (scaling) is usually observed. 
Despite the progress achieved in this field, the multi-species case is still poorly understood. Furthermore, when the 
H density of particles is high, the hard-core constraint on the dynamics of the diffusing particles becomes important. 
Experimentally the single-species fusion model used to describe the photoluminescence of excitons diffusing along one 
dimensional chains is seen to apply only as long as the initial exciton densities are small jjj. To our knowledge, there 
has been no systematic effort yet to investigate theoretically the regime of high densities of reactants. 
£j , This is the first in a series of technical papers where we will explore a new approach in an attempt to remedy the 
I* ■ difficulties encountered so far. We propose to start from a quantum spin chain formulation of the master equation, 
. £h ! fermionize and subsequently apply renormalization group techniques to deal with the interaction terms arising in the 
multi-species case. Note that in the single species case, the method applies at arbitrary densities as the hard-core of 
the classical particles is automatically accounted for by the Fermi statistics. When various species react and diffuse, it 
is appropriate to distinguish two cases: i) the various species have infinite on-site repulsion with themselves only; this 
can be treated readily following the methods outlined in the following sections; ii) the particles of different species all 
have a hard core constraint; this is far more difficult and will be investigated elsewhere. As with the other methods 
devised so far, there is a price to pay: the calculations involved are sometimes extremely tedious. 
The purpose of this first paper is modest as we will focus on the single-species case for which fermion-fermion 
interactions do not arise as a consequence of mapping classical particles to fermions. We want to show explicitly 
how by elementary means we can reproduce known results (the long-time behaviour of the density), and derive some 
new results, i.e., the explicit and exact analytical form of the two-point correlation function of one and two-time. 
The paper is intended as an introduction to the fcrmionic functional integration approach which we will combine in 
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forthcoming papers with the renormalization techniques to treat the multi-species case. 

In the first section, we define the model, introduce notations and review an extension of a method developed by 
Lushnikov. The second section is devoted to the application of an elegant technique introduced by Grynberg et al. 
in a different context to evaluate the two-point (one-time) correlation function. The third section deals with a new 
and powerful formulation of the problem in terms of path integrals of pseudo-fermionic variables. We illustrate the 
technique by computing the two-time correlation function. A brief discussion concludes the paper. 



II. LUSHNIKOV'S APPROACH 



This section introduces notations, summarizes and extends Lushnikov's genuine approach. We consider a lattice of 
N (even) sites (length L = Na, a = 1 and assume ^ even), with periodic boundary conditions, on which classical 
(spinless) particles with hard-core can diffuse (annihilate) to adjacent empty (occupied) sites with rate D. Whenever 
the arrival site is occupied, an annihilation-reaction (1 + 1 — ► 0) takes place. A source of intensity J injects pairs of 
particles on adjacent sites (0 — » 1 + 1). Lushnikov || has managed to rewrite the master equation that describes the 
annihilation and diffusion processes described above in terms of an imaginary-time Schrodingcr equation: 

j t \m)=mt)) , (i) 

where C denotes the so called "Liouvillian" , which by abuse of language we will call "non-hermitian Hamiltonian" . 
\il>(t)) represents the state of the system at time t, 



i^(*)>=E( p « n >,*) n i°) > ( 2 ) 

{"} m'({"}) 

where m'({n}) represents the sites of configuration {n} which are occupied. The Liouvillian is given by Q 

£ = (J + D)^2 (o'm^m+l + <T m°m+l + a m a m+l + ff m a m-l) + D X! (°m+l°m ~ a m+l a m ~ ^ a m a m) ~ JN ( 3 ) 



When J = 0, i.e., no source, in addition to diffusive processes with rate D, only irreversible reactions (1 + 1 — > 0), 
with rate 2D take place. 

For finite (J > 0) source, the diffusive processes (1 + 0^0 + 1 and + 1^1 + 0) take place with rate J + D and 
we also have reversible reactions: particles are annihilated (1 + 1 — > 0) with rate J + 2D and created (0 — ► 1 + 1) 
with rate J. It is worth emphasizing that these rates are not independent, and are chosen such that the Liouvillian 
is quadratic in the spin variables for a single species (the higher order terms cancel due to the properties of Pauli 
matrices). We also point out that in this model the "annihilation rate" (J + 2D) is always bigger to the "creation 
rate" (J). In the n-species case, this property does no longer hold if we assume hard core repulsion between all species. 
Indeed, one obtains a spin Hamiltonian (S = n/2) that is a polynomial of higher order in the spin operators and in 
general cannot be solved exactly. If however, we assume infinite on-site repulsion only bewteen particles of the same 
species, we can rewrite the Hamiltonian as a quadratic form of coupled spins 1/2. The latter can be solved by the 
techniques presented below. 

To solve the Schrodinger equation in imaginary time, Lushnikov performs a Jordan- Wigner transformation and intro- 
duces the fermionic operators a m = Ylj< m (^ ~ ^ n j) a rn- 

Because of the form of the resulting non-hermitian Hamiltonian, it is appropriate to work with Fourier modes 
a q = 7^ Em a m e~ lqm . The antiperiodic boundary conditions Q lead to: q = ±(21 - l)v/N , 1 = 1,2,..., N/2. On 
Fourier transformation, the evolution operator reads C = J2 q >o wnere ( n q = a ^ a q) 

C q = 2(J + £>){ cosq(n q + n- q ) + sin q(a 9 a_ g - a+a± 9 ) j + 2L> j smq(a q a- q + a+al g ) - (n q + n- q )\ - JN , (4) 



Now, by a BCS-like Ansatz, 



<J>0 ? >o 
Lushnikov || is able to decouple the dynamical equation as 



w)) = ni^w) = n(^(*)«+^w)io), (5) 
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j t \%(t))=£ q \%(t)) ,Vq>0 (6) 
For a lattice which is initially completely occupied, i.e., A q (Q) — l,B q (0) — 0, one solves the above equations by 

4( J + 2D) sin (|) 2( J + 2L>) sin (|) 

where 

pi = -2(J + 2£>)(1 - cosg) , p 2 = 2 J(l + cos g) (8) 
In the absence of source (J = 0), the solution simplifies considerably to 

A q (t) = exp(-4Dt(l - cosg)) , fl g (t) = cot (|)(exp(-4D*(l - cosg) - 1)J (9) 

At this point it is worth noting that the ket \ip(t)) characterizes the state of the system at any time without however 
being an eigenvector of C (The Liouvillian is not normal, however see below) . 

In ref. [Q], Lushnikov calculates the density of particles by the method of the generating function, which we extend in 
order to evaluate the two-point correlation function. The density reads 

p(i)=5Z^(W)P({n},t),Vi, (10) 
M 

Similarly, the two-point correlation function is written as 

g r (t) = {n t n l+r )(t) = J2ni{{n})n i+r ({n}) P({n},t) , (11) 
M 

where the translational symmetry of the system has been used. We observe that in this formalism 

p(t) = (0\exp(£<r-)ni\il>(t)) , Q r {t) = (0| exp(E a~) n i n i+f .^(t)) (12) 

n n 

as one can check using the explicit form of \i()(t)) (develop the exponential, order each term and perform a Jordan- 
Wigner transformation ||). It is appropriate to consider the following generating function: 

G(x,y,z,t) = (Q\exp(x ( rr + ya7 +r + z E a~) |V(t)) - E x^y^z^ ^-"^({n}, t) (13) 

n^i,i+r {n} 

So we have, 
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dx 



d 
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p{t) = —G(x,z,z,t) , Q r (t) = irir G{x,y,z,t) (14) 



x,y,z—\ dxdy 



x 7 y,z—l 



To compute the generating function, we rewrite the state as 



\m) = n { a s)« + B q (tj) io) = n (B q (t) e « sin ?(" - i°) ( i5 ) 



Note the normalisation condition due to the conservation of probability: 



e *) = n w - a « (*) ° ot |) = 1 

{n} g>0 



Next expand the argument of the exponential as 



G(x, y, z, t) = (0| ( 1 + (xai + ya,i +r + z E a «) + { z2 E a « a «' + xz(a.i E a « + E a " a ^ 
+yz( E a n a i+T . + a i+r E «») + zya i+r ai} + . . . J ]~[ ( 2 — a+a+ sing 

n>i+r n<i+r,n^i / q>0 \ n>m 



(n-m)-fl,(f) ]|0) (17) 
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In this expression, only the terms proportional to "xy" contribute to Q r (t). Let us call the first of these terms G\, 

= *4 ( £ cot r% ) ( n - ^(t) ) ( 18 ) 

In the absence of source, we have A q (t) — ► and B g (£) — > — cot | exponentially fast (see ^), so that in the thermody- 
namic limit (N — > oo), the asymptotic behaviour of (? r (i) follows as 

We can do the same for the density and in the thermodynamic limit, one obtains 

p(t) = JJ{0| ( 1 + (ma, + z^] a„^ + |z 2 ^ a„a„/ + aw; (a* ^ a„ + ^ a n ai) + . . . | J x 

g>0 \ n^i i^n.n' ,n>n' n<i n>i / 

x f E "ktf sin ?(" - m) ~ B q (t)) |0>| K _ f - l (t) ( 2 0) 

\ ™>™ / ' 01 A q (t) cot | 

In the above, we used the following identities: 

1 ^ ^ g i cosg(z - ^ . ^ \_^( J cosg(i — |) 



E ro >iSm ? (m-«) = -(^cot - + ^ ) , 2^ sin 9( m " i ) = 2V COt 



- - - sin § 



n>m 

The evaluation of the two-point correlation function requires the calculation of all the terms proportional to xy, which 
in general is a very hard task. In the following sections we will be able to solve this difficulty by reformulating the 
problem in a different language. 

Using the explicit form of A q (t) and B q (t) (0^) and the results of appendix A, we find the asymptotic behaviour of 
the density in the (irreversible) critical case as Q 

p -ADt rw i 

p(t) = / dqe* Dtcos « = e- iDt I (ADt) ~ -== (22) 

7T Jo VSnDt 

Similarly the asymptotic behaviour of the two-point correlation function is 

7T Jo sing 



= e- 4Dt 

0<n<r 

which implies 



E \hn-r+l(4Dt) - hn-r^Dt)\ ~ ^ (23) 



C ^ ~ -^m (24) 



Unfortunately, in the massive case (when the source intensity is finite), this method applies only to the computation 
of the density. Assuming Dt, Jt ^> 1, we find 

/oo 
dt , e -4(j+D)t' (j ( 4jDt ') _ i^ADt 1 )) (25) 

'i + ^r^™ (26) 



+ VJ + 2L> V 2£»/ 8JW8irDt 

:ium value of the density, in agn 
the next section we provide the full two-point correlation function when J > 



where p eq — -jj^^== represents the equilibrium value of the density, in agreement with Lushnikov's result . In 
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III. THE PSEUDO-FERMIONIC APPROACH 



In this section we evaluate the full two-point correlation function in the general case by means of a powerful 
formalism. The central idea is to perform on the fermionic non-hermitian (and non-normal) Hamiltonian a generalized 
Bogoliubov transformation which allows us to work with a diagonal evolution operator (see Ref. Jl^,|lJ]). Following 
previous works ]l3|-|l5|,|7|] , we denote each of the 2 N possible configurations by a ket |n): 

(n|n') = ««,»' . £l»X»l = 1 ( 27 ) 

n 

In this Fock space, we can efficiently record the probabilities for the various configurations in the ket 

\P(t)) = Y,P(n,t)\n) (28) 

n 

The master equation governing the dynamics of the annihilation and diffusion processes described in the previous 
section can be rewritten as 

J^I-PW) =W|P(t)) =5^ftP(n,t)|n) =J2( A (n'^n)P{n',t)-A(n^n')P(n,t)^\n) (29) 

n n,n' 

where U denotes the evolution operator, A(n' — > n) and A(n — > n') represent the transition rates J and _D in 
Lushnikov's formulation. The matrix elements for the operator are 

(ri\U\ri) = A(n -> rc'),Vra' 7^ n (30) 
n — * n') (31) 



This Fock-space formulation was used in Ref. |L3 14 to study a stochastic adsorption-desorption problem. In what 
follows, we will specifically focus on the reaction-diffusion problem in one dimension. 
Let us now introduce the left and right steady-states, respectively 

(X|=X>| , \x) = ^P(n,eq)\n) (32) 

n n 

where P(n, eq) denotes the probability for a configuration \n) at equilibrium. It is easy to check that e has no effect 

on |x) and (x|, and the conservation of probability leads to (xlx) = !• 

The transition probability from a configuration \n) to \n') is simply: W n>n '(t) = (n'\e \n). 

We intend to calculate the density and two-point one-time correlation functions of a system initially in the state 
\<j>a) = P(n,t = 0)|n). The occupation number operator n r being diagonal in the basis {|n)}, we have 

p(t) = J2(n'\nr\n'}W n , n ,(t)P(n,t = 0) = ^(n'Ke £ *|n}P(M = 0) = {x\n r e ct \4> ) (33) 

and similarly, 

Qr{t) = ^{n'\rnn m \n')W n , n ,(t)P(n,t = 0) = ^{n'\nin m e ct \n)P{n,t = 0) = (x\nin m e ct \cf> ), (34) 

where r = \m — l\. At this point, we perform a generalized Bogoliubov transformation (rotation supplemented by a 
rescaling), 



^ \ _ / a cos B q a 1 sin 9 q \ / a+ 
(,-q ) \ — asin6 q or 1 cos 6 



(35) 



in order to obtain a diagonal representation for the evolution operator. 

This transformation is orthogonal, i.e., invertible and canonical, in the sense that it preserves the anticommutation 
relations of the a q 's, namely: 

{# = = = (36) 
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Despite the fact that the £ g and are not adjoint of each other, this representation will be very useful in the following. 
We set a — ( j+ 2D J , so that the mode q evolution operator becomes 



-C q = 2[D(1 - cosg) - Jj [(^ q + £ q £_ q ) + S m6 q 2 (Z_ q £ q + + 2(£t ff # + 

+ y/J(J + ^D) sin<z(cos20,(£+,£+ + + sin2fl,(^+ + + 2 J (37) 

To get rid of the terms that do not conserve the number of pseudo-particles, we choose 9 q as 



tan20 9 = ^— £ 38 

(J + D) cos q — D 



so that the Hamilton operator becomes 



c = - E A ? + = - E & ( 39 ) 

9>0 q 

where 

A 9 = 2 (.0(1 -cosq) + j) (40) 



on account of the periodic boundary conditions (X) q >o cos< ? = Now it is readily seen that (\\ and \x) act, 
respectively, as left and right vacua, i.e., £ q \x) = and (x\£q = 0- To simplify the calculations, we express the initial 
ket-state \4>o) in terms of the steady state |x). We consider here two kinds of initial conditions: 

i) The whole lattice is initially filled. We write |</>q) = \all) and immediately conclude a q \all) = 0. Using the inverse 
of (||^), one can check that 

\all) = JJ (l - cot 8 q ^e- q ) W) = exp ( - £ ^%£± 8 ) Ix) (41) 

ii) The lattice is initially empty; one can check in the same way E3 that : 

l</>o> = |0> = (l + tan^+elj |X> = <*P (£ £,) lx> (42) 

g>0 g 

The time dependence of (i) and £^ (i) follow as 

&(*) = e" £t e 9 e £t = e- A ' 4 a (43) 
e fe + W = e-% + e £t = e A ^+ (44) 

(45) 

= 0) = (xl&iffalaiO = cot6 k J kl: - k2 (46) 



Futhermore, we have 



and 

(£,k 1 £k 2 £k 3 £,k i )(t = 0) = cote» fcl cot0 ka 5 kl ,-k 2 Sk 3 ,-ki + cot^ cot6»fc 2 (<5 fel ,-fc 4 4 2 ,-fe 3 _ 4 1: -fe 3 4 2 ,-fc 4 ) (47) 

as one can check by applying Wick's theorem (see also section IV). Using the properties of Fourier transform and 
the generalized Bogoliubov transformation ([35]) , the expression of the density and the two-point correlation function 
become respectively (for a lattice initially filled, |</>o) = \all)): 

p(t) = ^ E sin ' 6k - E s[ndk cose kl (x\^-k^e Ct \all) (48) 

k k,k' 



G 



and one can do the same for unconnected one-time correlation function Q r (t) (|34|). 

To derive tractable formulas, we have performed tedious but straightforward calculations. Indeed, we have extracted 
the time-dependence of p(t) and Q r (t) using ( p3| ) and commuted all the pseudo-creation operators to the left of the 
pseudo-annihilation operators. In the expression pit) and Q r (t) only terms like k£k' MO — (£— feCfc')ft = 0) and 

3f (|| g). 



(xl£-fc£fc'£-q£g' MO — (£,-k£,k>£,-q£,q'}{t = 0) survive: these were evaluated with the help of 
In the thermodynamic limit, we arrive at 

p(t) = 22 ( sin 2 9 k + e~ 2Xkt sin 6> fe cos 9 k cot fc 



N 

k>0 



^ £ dq(sm 2 d q + e- 2 ^ 1 cos 2 9 q ) = ^_ + ^L_= + 2(J + 2£>) j™ dt'(l Q (ADt') - h{ADt'))e- il - J + D ^ (49) 



which coincides with (|2E 

It is worth emphasizing that this result is general and works for both, the massive (J 7^ 0) and the critical (J = 0) 
cases. The point here is that the limit J — + is not singular, despite the divergence of cot 6 — > — 00. In fact, integration 
over k and k' of terms proportional to sin(9fc cot 9k> yields finite results. Therefore, we can perform the computations 
([35]) at J finite and set subsequently J = in p(t) and Q r {t)- 

Similarly, the two-point correlation function (of one-time) is evaluated as (r = \m — l\): 

Gr(eq) = pl q -\ j \J dq sin 2 9 q cos qrjy J dq' cos 2 Q q ' cos q'rj + J dq sin 29 q sin qr J (50) 

</r(*) - £r(e<2) = (p(t) 2 - p 2 g ) + ^2 ( _/ ^ e COs2 « C0S ( / COs2 6 l' C0S 9 ' r 

+ — ^ ( y ^9 sm sin qr) ^ y dq' e~ 2X i' 1 cos 2 6>q' cot 8 q > sin q'r 

e -2A,j,t cog 2 cosq ' r 



i ^ y dq sin 2 6* g cos qr) ^ ^ dq' 
^~ 2^ y dq sin 29 q sin qr*J J dq' e~ 2 Xq,t sin 29 q > sin q'r 



— / dq g~ 2 ^<j* cos 2 0„ cos qr ] - ( — x / dq e~ 2A<jt cos 2 0„ cot 0„ sin qr 
Wo / Vtt 2 7 

dq e~ 2 Aljt sin 2 6» 9 cot 9 q sin qr) , (51) 



where we separated the static contribution to the correlation function from the dynamic one. Using known properties of 
modified Bessel functions given (see appendix) and writing I' m (x) = -^I m {x), we finally infer (notice that Q r (eq) = p 2 q ) 



G r {t) = p{t) 2 + 2(J + 2D)(A - Co) ^ dt'e- i{J+D ^[l r {ADt') - I' r {ADt')} 



^/J(J + 2D)B Q 



dt'e-^+W^—UADt') 



2 J t 2Dt 

J + 2D 



*J I! , 



poo 

£ ( / dt ' ^ i(J+D)t {2J 2n _ r+1 (4Dt') - 2/ 2n _ r (4Di') + l' 2n _ T {ADt')} 

<n<r Jt 

poo 2 

(2(J + 2Z?) / dt' e -^ J+D ^' {l r {ADt') - r r (4Dt')}) 
2(J + 2D) 2 ( ^ ^-r)(rdt'e-^^' l2n -^ Dt ' 



0<n<r 



J2 (y°° dt' e-^ J+D »' {2I 2n _ r+1 {ADt') - 2I 2n _ r {ADt') + I' 2n _ r {ADt')})) (52) 



0<n<r 

In the above formula Aq 7 Bq, Cq have been defined by 
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If 7 ' If 7 ' If 7 ' 

Aq = — dq cos 2 9q cos qr, Bq = — I dq sin 2(9 g sinqr, C$= — I dq sin 2 9 q cos qr (53) 

77 Jo 77 Jo ' 71 Jo 

or more explicitly, in the massive case (with help of (A12)), 

A „ . i£±ia r f (c089r _ e^l±i) ^.(r-D ) . _ {J + 2W .-, o-gL, (54) 

7T Jo A, V 2 / 4V J 2 + 2 JD 



Bo = jT g ( oo. g (r + 1) - cos 9 (, - 1)) = -^^A (55) 

7T Jo A, V 2 J 4V J 2 + 2J/J V ' 

with the notation 

C = — < 1 (57) 

(J + D) + V J 2 + 2 JD " 

Note that the above result, though obtained by elementary means, is new. The formula can be evaluated numerically. 
The connected two-point correlation function C r (t) = Q r (t) — p(t) 2 follows immediately from the above. To establish 
a connection with the results derived previously w\ we explicitly evaluate C r (t) when J — 0. Taking due attention to 
the apparent singularities occuring in this limit, we find 

1 f 7 ' C 
C r (t) — — ~ / dq sin 20 q sin qr / dq'e^ 2X i' t cot9 q > sinq'r 
27r 2 Jo Jo 

' 1 f 7 ' \ 2 / 1 f n \ / 1 f" \ 

-J dq e ~ 2Xqt cos qr) ~ \~ J d 1 e~ 2Xqt cot 6 q singrj (- y e ~ 2A<!t sin 2 cot 0, sin grj (58) 

or in terms of elementary functions, 

£.(*)= ^ e- 4Dt (/ 2n _ r+1 (4/Jf)-/ 2n _ r (4/Jf) 

0<n<r 

+ ( £ e- 4 ^/ 2n _ J ,) 2 -( £ e-^I 2n _ r+1 ) 2 -(e-^I r {ADt)f (59) 

0<n<r 0<n<r 

This expression is equivalent to the result derived by a well-known mapping of Lushnikov's model to Glauber's ID 
Ising model 0. 

Finally, we compute (using formula (All)), the asymptotic behaviour of the two-point correlation function C r (t) when 
Dt > 1, Jt > 1, and r < oo 

C r (t) ~ _ C r - X e^ ((j + 2D){1 - C) 2 + J(l + C) 2 + rJ ( J + 2J> ) (i _ C 2 )| - 1 = (60) 
8Jt\/8^Dt\ -D M VJ(J + 2Z?) 



As expected for a finite source intensity J > 0, the density and one-time correlation function decay exponentially 
with time. Morever, the one-time correlation function also decays exponentially with distance. 

A similar calculation can be performed for the case of an initial empty lattice ( 1 Q ) = |0)). In this case the density 

is p(t) = 2J J t °° dt' (7o(4£>i') + h{ADt'))e-^ J + D ^' ~ p eq - e - 4Jt /V8^Dt @. It has been shown @, via the 

mapping to the Glauber dynamics, that the nearest neighbor connected function decays (Dt, Jt 3> 1) as C(r — 1, t) ~ 
p eq e- 4Jt /V2^Dl. 

The large-time and large-distance behavior of the one-time correlation function can also be studied. In fact, in the 
critical case it has been shown [^|J^] that the one-time correlation function obeys a scaling form. 

Here we investigate the behavior of C r (t) in the limit where r and Dt — > oo with u = finite. In the massive case, 
we do not expect a scaling form : 



Cr(t)~-\ll + 2 4C-\l-< 2 ) 9 -^fi, (61) 



J v s ' Ar 2 

where 
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The effect of the source is to disrupt the spatial correlations, i.e. to make them short-range. In this sense the finite 
source prohibits the formation of arbitrarly large vacancy domains. 

In the next section we will introduce a field-theoretical approach to deal with two-time correlation function. Our 
approach and the results obtained thereby complement the exact treatments of the critical case jig] . It is worth 
pointing out, that Glauber calculated in his pioneering work the two-time two-points of spin correlations functions 
II ■ 



IV. FIELD THEORETICAL REFORMULATION 



The purpose of this section is to reformulate the results of the previous section in a field theoretical language, 
i.e., in terms of path integrals of fermionic variables (see for example fin\). We define the Grassmann numbers 
rjq,rf" q which anticommute with each other and with the fermionic operators introduced in the previous section, i.e., 

{Vq,Vq'} = {Vq,Vq>} = {Vq,Vq'} = and iVg^q'} = {Vq, Cj} = {Vq,$q'} = {v* q ^q>} = 0. We follow standard practice 
and consider the coherent states associated to the fermionic variables. We recall that pseudo-fermionic operators 
required the introduction of a left vacuum (x\ and a right vacuum \x), ergo we define the right \r]) and the left 

coherent states (% respectively, as = e q q \x) , (v\ — e q 9 9 

Despite the fact that £ and £ + are not adjoint of each other, we find the familiar results: 

m=W = l , m=<£S r >* and £ q \r ] )=r)q\v) , fflf, = ^{v\ , = {rj\t g , ^\v)=-^\v)- 

Most importantly, the closure relation holds true, i.e., 

l[dr ] * q dr 1q e-^ r, >«\r)){rf\ = t (63) 
q 

At this point, we know from field theory how to calculate (x\n m e ct \all) and (x\n m n m+r e ct \all) using the path-integral 
formalism. We discretize time in M infinitesimal intervals of width e = limjv/^oo jj . C is normal ordered, and the 
closure relation is inserted into the above formulas. We have for the density 

p(t)= f [] dr,l a dr, q , a e-^<^{x\n m \mi)^ , (64) 

«* „ 1 n ji 



q.a=l...M 



where 



lim (^|e £ei |^0 = e ^ q cote i<^t- q ^ 

M -^oo 

= ^(X!™ 2 ° k ~ ^2$' m 8k COS 8 k >(x\£-k£k>\llm)^ 
k k,k' 

jyfy^ sin2 ® k ~ E Sin ® k C0S ® k ' >■'!>■! >■) ( 66 ) 



N . 



Taking the continuum limit ]l7|] , we arrive at 
p(t) = ^J n^W^W(E sin2 ^^E sin ^ cos ^' 1,(t)v4" E '^'' (0l ''' !0) e" So!,,(t) '' (t)) , (67) 



k.k' 



where, the Liouvillian operator C acts as the Hamiltonian of Field theory. By abuse of language we call So(n* (t), rj(t)) 
the Euclidian action of our problem, which is defined by 



S (v*(t),r,(t))= / dt' E(<( f ') 5 t'^(«')-A^(i'),^(O))+E^'(°)v(0) (68) 

J ° q q> 

Taking averages (. . .)s on the gaussian distribution, the density can be rewritten as 







p® = jf ( ( E sin2 0fc - E sin * cos % w) e ~ E < ^ (0) "" 9 (0) > 

We proceed to discretize the free Euclidian action So which is bilinear, as 



So 



/ 1 0\ 

-a 1 

-a 1 ... 



S(q) 





V 



-a 1 / 



where, a = 1 — -o> A 



Following Ref. (1J, we find det(5(q)) = 1 and 

lim S-] 3 {q) = e-^t*-"*"), Va = 1 . . . M > /3, t c 



at 



(69) 



(70) 



(71) 



Futhermore we have, in the continuum limit ^ / IIg d"rj*(t)dr) q (i) e s ° = lj 

(»fe(fe)»£(*i)>Sa = / n<W) d »fc(<) »b(*a)^(*i)e- So = We-^^ - * 15 
Applying Wick's theorem, leads to 

<^> (*) = (%(i)V 1 E * COt M « ) c = S p ,. p , cot ( 

* ' On 



(72) 



(73) 



With this result, it is straightforward to calculate the density using (|73| ) and (|6£ 

The computation of correlation function requires evaluation of terms like (£/ci£fc2£fc3£fc4) (t) . We generate in the 
standard fashion the multi-point correlation functions. The generating functional in discretized time reads: 



Z \9q,a', 9q' a '] 



J2„ {~5 cote vVliV- q ,i+Y /a > 2 ( -9" q ,<>V q , a +v" q , a g q ,c.)} 



(74) 



S 



here g q ,a,g qa denote (Grassman- numbers) the coefficient of the source terms. Note that this functional differs from 
the usual field theoretical one ]l7| by a term that codes the initial condition, i.e., Y[q>o eC ° teqVq ' lV ~ q ' 1 ■ Taking this 



term into account however is no trouble, i.e. 



M 



zig q , a ;9* q ,, a/ }=\ n( e9i - is;icote5 +3- 9 .i.9;i3-,,i5 9 ,i) n n e^^^V.A 

I q>0 ) q,q' i,j-2 



(75) 



Considering that the source term will be set to zero in the end of the calculation and noting that for a = 1, only 
pseudo-creation operators contribute, we find it convenient to work with 



q 

From the definition, it follows then that 
or more generally, 



(76) 



S 9-qA S 9* q ,l S 9* kl ,M S 9*k 2 ,M S 9q,2Sg-q,2 



(77) 



g=g'=o 
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5^Z[g q , a -g* ql ^\ 



^*,„,iK,i • • • S 9- qi ,i s 9* qi Jg* kuM Sg* k , iM . . . 5g* knM 5g* KM 5g qu2 5g- qi .2 . . . 5g qn , 2 dg*_ qnfi 



(78) 

3=9* =0 



As an application of this field-theoretical formulation we will now extend our approach to the computation of the 
two-time correlation function, for an initial state \(f>o), 

G r {t,t') = (n J+r (t + *>,(*)) = (x\nie ct ' n l+r e ct \M = ^(n> ;+r e £ *'n ; e £t |n)P(n, 0) 

n,n' 

= (n'\n l+r \n')(n''\n l \n'')W n/n 4t')W nl , n (t)P(n,0) (79) 

n.n' ,n" 

As above, we consider an initially filled lattice (\<fro) — \all)). 

To perform the computation of Q(r, t, t') we need to evaluate terms like 

(xl^'^^VO = ^ (80) 

This expression can be calculated by using the generating functional, where as before t (t r ) is discretized into M (N) 
infinitesimal time steps (respectively), with M, N — > oo: 



l~\C e Ct' c <■ Cti ;;\ S 12Z \9q,a\g* q ' a >] 

(xKpCp'e £ q £q>e \all) 



S 9% 2 Jg^ A Sgl qiA Sg* qul Sgl M Sg* qliM Sg; M+N 5gp M+N Sg qu2 Sg- q ^ 
e -(^+V)(*+*')-(A,+A,,)t ( cot 6p cot e q s p ,- p ,s q ,- q , 



g=g'=0 



+ COt 6 p COt 6 p > (j)-p^S- q ^ - 5- qt pS- p ^ q '^j 



(81) 



where the continuum limit for the time has been taken. 
Applying the same technique to each term of Q(r, t, t'), yields: 



G r (t, t') = p\ q + 2(J + 2D) Peq J™ dTe-^ J+D ^ T (l Q {ADT) - h{ADT)) 
+ J(J + 2D)^£ dTe- 2{J+D ^ T {j r (2DT) + I' r (2DT)^j 

*(f dT ^ +D)T {^-)-^(-)}) 2 
+ J{J \ 2D) ^fdTe-^ T ^I r (2DT)J 

- 2(J + 2Df(^J™ dTe- 2(J+D)T {l r {2DT) - I' r (2DT)^j 

x ^jf dTe-^ J+D ^ 2T+t '^I r (2D(2T + t')) - I' r (2D(2T + t'))^ 

x ( E /"^ e - 2 ( J+D )( 2T+t '){2/ 2 „_ r+1 (2^ 

\.0<n<r"'' / 
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X 



X 



J( J + 2D) ( J™ d Te- 2{J+D){ - 2T+t,) (l {2D{2T + <')) - h(2D(2T + i'))) 
^jf dTe- 2 ( J+jD ) T (/ (2Lir) + / 1 (2Lir))^ 

J( J + 2D) ^ j~ rfTe _ 2(J+0)(2r+0 _ fr(2 J(2r + ^ ^» dTe^ J +°> r ^/ r (2JXT)) 

2J(J + 2D)(j^ dTe- 2{J+D)T {l r (2DT) + I' r (2DT)^j 

dTe -2(J+D)(2T+t') |j r ( 2 £)(2r + i')) - I' r {2D(2T + t'))^j 
4( J + 2D) 2 f jf dTe -2(J+D)(2T+t>) (7 ( 2jD (2T + £')) - h(2D(2T + t'))) 
2(J + 2L>) 2 ^ jf dTe- 2 (' /+I5 )( 2T+ *'){/ r (2L>(2T + if)) - I' r (2D(2T + *'))}J 
2(J + 2D) 2 ( £ (2n-r) / 



2(J+D)(2T+t') hn-r{2D(2T + t )) 



e -2(J+-D)(2T+t') 



X 



D(2T + f) / 

/DO 
dTe -2(j+DK2T+t ) {2l2n _ r+1 (2D(2T + f) - 2I 2n _ r {2D(2T + t')) + I 2n _ r (2D(2T + t'))}) (82) 



0<n<r 

In the critical case some care is required when taking J — > 



r (t, = e - 4I? ( 2t+t ') I (j„(20(2t + If))) 2 - (l r {2D(2t + t'))) J 

+ £, e -2£»(2t+t')| J™ dTe -2DTj_ Ir{2DT) ^ 

x | ^ (j2 n _ r+1 (2Z)(2i + 0) - hn-r{2D(2t + i'))) I 

+ e" 4I W)j( hn-r{2D(2t + t'))) 2 - ( ^ 7 2 „_ r+1 (2D(2t + i'))) 2 } (83) 

I 0<n<r 0<n<r J 

To study the behavior at finite distance of the two-time correlation function C r (t,t') = G r {t,if) — p(t)p(t') , we 
distinguish the massive and th e crit ical regime. We begin with the massive case and consider Dt, Dt' ^> 1, Jt, Jif ^> 1 



and r < oo. With help of ( A8- All), we have 



J 



Gr(t,t')~p 2 eq +{l + —)p eq 
J \ e 



2DJ r q 8JtV8^Dt 
i , p-iJt' 

C r (i,i')~-(1 + — )Pe q ; (84) 

rv I y 2DJ Feq 8Jt\/8^DF V ' 



8Jt'V8^Dt' 

When r < oo, Dt and Jt are finite and Dt' , Jt' > 1 we obtain: 

J + 2De- 4Jt ' J + 2D e -^(*+*') 



Gr{t,tf) pi + 



128J (Dt') 2 4JZ? 16tt j DV< 7 (2< + ^) 1 

C P (*,t , )~-(l + ^K)Pe« — ; (85) 
rv j y 2DJ Heq 8 Jt'V8^DF V ' 
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In the critical case (r < oo), both asymptotic behavior Dt, Dt' ^> 1 and Dt finite with Dt' ^ 1 of the disconnected 
and connected correlation function are given by: 



r / 1 / 2t 



- 2 ■ < ? 2t\ 

4TrD 2 (2t + t') 2 V ' 8 V ^ ' T) 

C r (t,t') 3_= (86) 

8nDVtF 

Near the initial state (the density of particles is high), when Dt, Dt' <C 1 the decay is linear and independ of r, i.e., 

g r (t,f) ~ l-4D(2t + t') 

C r (t,t') ADt (87) 

We now provide a scaling form for the two-time correlation function C r (t,t'). It is known, from the duality with 
Glauber's model, that the single-time correlation function obeys a scaling form for large-time and long-distances 
{Dt and r — > oo with the ratio r 2 /Dt held finite The scaling form is found as C r (t) ~ r~ 2 f (r 2 / ADt) , where 

the exponent —2 is believed to be universal JtJ. We further assume the long-time and large-distance scaling form 
C r {t,t') ~ r~ v h(u,v), where r, Dt, Dt' — > oo with u = r 2 /ADt and v = r 2 /ADt' held finite and arrive at 



C r (t,t') 



^i.Kh- K (ynerty/v-2VK) +K(l 



-IK 



where < K = < oo. 

At this point, it is appropriate to review what we have achieved so far: we have been able to reformulate the problem of 
the evaluation of the multi-point correlation functions in a language that parallels the field theoretical one. This allows 
us to compute in an efficient and systematic way physical quantities of interest despite some technical differences to 
the standard approach. While this paper deals with a free "field theory" of pseudo-fermions, it is tempting to apply 
the same formalism to the multi-species case where two- or multi-body interactions arise. The latter however will be 
investigated in a future work by pertubative renormalization group techniques, as no exact solutions are available. 



V. CONCLUDING REMARKS 



We have studied three different approaches to the problem of diffusion-annihilation of classical hard-core particles 
moving on a one-dimensional ring. Though Lushnikov's contributions to the problem are genuine and undisputable, 
we have shown how an extension of his generating function method to evaluate the two-point correlation function can 
be cumbersome in practice, even in the simplest case available of a single-species. We have seen that it is advantageous 
to apply a generalized Bogoliubov transformation first used by Grynberg et al. |l3|Jl4|| in a different context. The 
evolution operator can be diagonalized, i.e., expressed as a quadratic form of two operators that are not adjoint of 
each other. Despite this fact the formalism resembles the standard one and appears as a powerful tool. Indeed, 
we were able to compute for the first time the full one-time and two-time correlation functions for an initially fully 
occupied lattice (other initial conditions can also be studied) in the presence of a finite source. We derived a scaling 
form for the two-time correlation function. We used the results of the first section (algebraic decay in the critical case 
and exponential in the presence of source) to check the asymptotic behavior of the density and two-point correlation 
function. We discovered that while in the absence of source, the modes at long wave lengths fully control the long-time 
asymptotics of the density and correlation functions, in the presence of a finite source all modes contribute. This 
means that in the general case, the long wave-length approximations that were so successfull in strongly correlated 
systems (such as bosonization or conformal field theory techniques) do not work for the problem at hand. Moreover, 
the idea of exploiting the integrability of some spin Hamiltonians on which the multi-species Liouvillian maps, might 
turn out more elusive than expected. In view of the above remark, it would seem extremely difficult if not impossible 
to extract from the exact Bethe-ansatz solution of the non-hcrmitian spin Hamiltonian the relevant matrix elements 
that in turn allow the evaluation of correlation functions |7J]. We propose to tackle the multi-species problem in 
terms of fermion functionals, the main difficulties arising from the two- and/or multi-body interactions occuring in 
the process of mapping classical particles to fermions. We will illustrate the power of the formalism in a forthcoming 
paper, where we will apply the renormalization group scheme. Besides the obvious advantage of formulating the 
problem in a field theoretical language (perturbation theory,...), the method is applicable to arbitrary densities of 
particles, and thus complements the approach developed by Cardy and collaborators p~pr[. In higher dimensions, 
however, Fermi statistics requires the introduction of a gauge field that is strongly coupled to the fermions. We also 
intend to explore this line of research in the future. 
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APPENDIX A: USEFUL RESULTS 

In this appendix we provide some useful properties of the Bessel functions of imaginary argument. We recall the 
definition of the modified Bessel I n {z) function (n integer) 

(—i) n f n 1 r" 

Uz) = ± — '— e zsin{p)+lnp dp= - e zcosf> cos {n6)d6 (Al) 

2?T J_„ 7T J 



with I n (z) = I- n (z). We also use in ( [52| , [S2| , j83| ) the well-known properties I n -\{z) — I n +i(z) — ^f-I n (z) and I n _i(z) + 
I n +i(z) — 24zl n {z). The following integrals occur in the evaluation of the two-point correlation function. 

J x (r,t)= r^^e 4Dtco ^dq , h(r,t)= f ^ cos q e 4Dt cos Uq , I s (r,t)= f ^ cos2g e 4Dtcos ^dq (A2) 
J sing J sing J sing 

Setting q = q — it, e real and e > 0, we have 

smgr =Um singr 
sing ej.o sing 

^SL = lim Hm y( e -^(n-(^))_ e -2Mn + (^))) (A4) 

sing el ~^ q sing e|o,g-^^ ' 



Therefore, 



I 1 {r,t)=r( lim Vfe-^-^-e- 2 ^^))) ] e 4C *- s «dg 



/0 V e 4-°>9— ?„>() 

= tt ^ / 2n _ r+1 (4£>t) (A5) 



0<n<r 



Similarly, we find 



J 2(M) = | J! {-f2n-r(4I>t) + /2n-r+2(4Dt)}=7r hn-r^Dt) (A6) 



< n < r 0<n<r 



^M) = f H {h n -r-\{±Dt) + I 2 „+r-l(4Dt) + hn-r+^Dt) + I 2n +r+3 (4Z>t) } = 

= tt ^ 7 2 „_ r _ 1 (4Dt) (A7) 



0<n<r 



0<n<r 



With help of asymptotic behaviour of Bessel functions and using the properties of the Incomplete Gamma functions 
ijisi , the asymptotic behaviour (Dt ^> 1, Jt ^> 1) of the following integrals is readily found 



/ dt' , ( l - -L + 0((Jt)- 2 ) ) (A8) 

Jt V8^D¥ uVS^Dty 8Jt U ' ') K ' 



dt '7m^~iT^(m-sh;. +0((Jtr3) ) (A9 » 
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oo p-iJt' p-iJt . 1 v 

dt' ( -— 7 + 0(( Jt)- 3 ) ) (A10) 

A futher result (Dt 3> 1, Ji ^> 1) used in the evaluation of the asymptotic behaviour of the density and the two-point 
corr elation function (||,^,|l],||,|8|,||,||,|86|) 

-4,/t 



r°° , P -4jt , i 



+ (" 2 -i'<" 2 -3)(l28W + « Jt )" 3 ») + -} < A "» 
Finally, the calculation of the constants Ao, Bq, Co, when J > is performed with help of the formula Jig] 

^ 

V^a 2 -p(a + A/a 2 - /3 2 ) 



e- ax I u {J3x) = —^=^77, WRe(u) > -1, i?e(a) > i?e|/3| (A12) 
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